First principles derivation of NLS equation for BEC with cubic and quintic nonlinearities at 

non zero temperature. Dispersion of linear waves. 
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In this work we presented a derivation of the quantum hydrodynamic equations for neutral bosons. We 
considered short range interaction between particles. This interaction consist binary interaction (/(r^, r^) 
and three particle interaction f/(ri,rj, r^), the last one does not include binary interaction between par- 
ticles. From the quantum hydrodynamic (QHD) equations for Bose-Einstein condensate we derive non- 
linear Schrodinger equation. This equation includes the nonlinearities of third and fifth degree. It is at 
zero temperature. Explicit form of the constant of three-particle interaction was taken. First of all, de- 
veloped method we used for studying of dispersion of linear waves. Dispersion characteristics of linear 
waves were compared for the cases. It were of two-particle interaction in approximation third order to 
interaction radius (TOIR) and three-particle interaction, at zero temperature. We consider influence of 
temperature on dispersion of elementary excitations. For this aim we derive a system of QHD equations 
at non-zero temperature. Obtained system of equation is an analog of well-known two-fluid hydrodynam- 
ics. Moreover, it is generalization of two-fluid hydrodynamics equations due to three-particle interaction. 
Evident expressions of the velocities of the first and second sound via the concentrations of superfluid and 
noncondesate components is calculated. 



I. INTRODUCTION 

At the theoretical investigation of Bose-Einstein conden- 
sation (BEC) with account of the three -particle interac- 
tion (TPI) the nonlinear Schrodinger equation (NLSE) is 
used. This equation contain cubic and quintic nonlineari- 
ties |liL2]- For describing inelastic scattering imaginary 
part of interaction constants is used in [3]. Three-particle 
interaction leads to forming of dimers. This process leads 
to loss of atoms from BEC 14,JJ. 

Opportunity of obtaining of the Gross-Pitaevkii (GP) 
equation from the microscopic many particle Schrodinger 
equation is demonstrated in |'^. The GP equation is the ex- 
ample of NLSE. In [7J the direct derivation of the GP equa- 
tion from many particle Schrodinger equation has done. 

In article authors derived NLSE for system of bosons 
being in the state of BEC, with the short-range interac- 
tion potential, in third order in interaction radius. If one 
keeps only the term of the first order in interaction radius 
then one gets well-known GP equation from many-particle 
Schrodinger equation. There taking into account only po- 
tential of binary interaction. At derivation authors do not 
use suggestion about density of considered system. There- 
fore, there occur interaction whether by means of scattering 
or during finite interval of time. Presented in 10] derivation 
does not exclude interaction of a few particles at the same 
moment. 

In this article we use the same method for studying of 
both a three particle interaction (TPI) and influence of tem- 
perature on dynamic of neutral Bose particles. 

Fundamental and detailed distribution of influence of 
non-zero temperature on BEC dynamic is presented in pa- 
per [8], where were considered the two-particle interaction 
only. The influence of the noncondensate atoms on BEC 
dynamic is considered in fS], but dynamic of noncondesed 



atoms do not included. Interference of superfluid and nor- 
mal components usually described by two-fluid hydrody- 
namic [9]- Present day examples of using of two-fluid hy- 
drodynamics can be found in papers [10-13]. Connection 
of two-fluid hydrodynamics with the kinetic equations is 
presented in papers In cited papers authors used 

Boltzmann like kinetic equation with the collision term. 

Due to studying of three particle interaction the ground 
state energy density in second order in \J pa^ is calculated 
in lfl5l - l2lll . where p is the particles density, a is the s-wave 
scattering length. The first calculation of the constant under 
the logarithm is presented in [15]. The constant depends 
on a and from the parameter described the lower energy 
3 — )• 3 scattering of the particles. 

The explicit three-body contact potential for a dilute con- 
densed Bose gas is derived from microscopic theory in 
I 22ll . The derivation is based on the quantum expectation 
values of products of single mode annihilation and cre- 
ation operators. The three-body coupling constant exhibits 
the general form predicted by Wu [19]. It depends on s- 
wave scattering length a, defined via binary potential like 
in [19]. For describing the properties of BEC at zero tem- 
perature, with three -body interaction, in [22] there were 
used nonlinear Schrodinger equation contained nonlinear- 
ity of fifth degree. 

Here we present some physical effects there is in BEC 
due to TPI. Stability properties of BEC with TPI are con- 



sidered in 



2511 ■ The modulation instability of the 
BEC trapped in an external parabolic potential, is investi- 
gated in ll23ll . The explicit time-dependent criterion for the 
modulation instability of the condensate was established. 

The influence of initial conditions on stability of BEC 
was studied in [3] using a Gaussian variational approach 
and numerical simulations, in three-dimensional trapped 
BEC. Abdullaev et.al. [3] discussed the validity of the 



2 



criterion of stability suggested by Vakhitov and Kolokolov. 
In the works iH |23l] for studying of properties of BEC 
with TPI there were used nonlinear equation with cubic and 
quintic nonlinearities. 

Abdullaev et.al. [24] discussed localized ground states 
of BEC in optical lattices with attractive and repulsive TPI. 
For this aim a quintic nonlinear Schrodinger equation has 
used. In ll24ll the existence of unstable localized excitations 
which are similar to Townes solitons of the cubic nonlinear 
Schrodinger equation in two dimensions is shown. 

Marklund et. al. |25] found the Vlasov-like equation for 
BEC. They used Wigner function and NLSE with cubic- 
quintic nonlinearities. 

Under the condition, that two-particle scattering length 
tends to zero, it may be realized by manipulating external 
magnetic field, far off-Feshbach resonance i26ll . T tends 
to zero too. In this situation main role has TOIR and TPI. 
In this case TOIR and TPI are compared. 

There are various generalizations of GP equation. Here 
we interesting in nonlocal NLSE. We present brief review 
of known nonlocal NLSE generalization of GP equation. 
For studying of BEC the NLSE with nonlinearities con- 
taining spatial derivatives of wave functions is used. The 
NLSE for BEC with nonlocal nonlinearity was obtained in 
ll27ll . In that article NLSE contains second spatial deriva- 
tives of square of wave function module. In |7] NLSE in 
third order on interaction radius (TOIR) approximation is 
derived. In this case NLSE is integro-differential equation 
with spatial derivatives of wave functions. Comparison be- 
tween nonlocal nonlinearity ll27ll and TOIR 0] for two 
dimensional space is considered in ll28ll . Local Lagrangian 
density, for two-body and three-body effective short-range 
interaction, contained dependence on (V'*^)^, and 
V{tlj*Tp)'S/{ip*iJj) is presented in ll29ll . There are given a 
determination of the strength of the three-body contact in- 
teraction for various model potentials. 

Generalization of GP equation to the large-gas- 
parameter regime is suggested in llsoll . It is NLSE with 
nonlinearities of fourth degree. 

In this article we use method of quantum hydrody- 
namic (OHD) of many particles system developed in works 
lE, 31, 22] ■ In llsill equation of QHD of charged particles 
system was found in external electromagnetic field taking 
into account Coulomb interaction of particles. In lf32h for 
system of charged and spinning particles there were de- 
rived equations of quantum hydrodynamic, i.e., continuity, 
balance of momentum, magnetic moment and energy. In 
isill there were developed methods of calculation of many- 
particle functions arising in equations of QHD. In |7] there 
was developed method of quantum hydrodynamics for sys- 
tem of bosons, fermions and mixtures, with short-range in- 
teraction. In this article there were given corresponding 
equations of continuity and balance of momentum. On the 
base of system QHD equations the NLSE describing the 
dynamics of boson-fermion mixtures is obtained. Particu- 
lar case NLSE for bosons noticed in |0] is the GP equation. 



For boson systems with two- and three -particle interac- 
tion, the problem of finding the quantum stress tensor can 
be solved and conditions of its existence can be clarified. 
Here, we will limit ourselves by the development of the 
quantum stress tensor and show that under the standard as- 
sumptions this tensor can be transformed so that the mo- 
mentum balance equations for bosons coincide with the 
analogous NLSE contained cubic and quintic nonlineari- 
ties. Using this momentum balance equation derivation, 
we can obtain the equation (for wave function in medium 
or the order parameter), which is NLSE with nonlinearities 
of fifth degrees. Quantum stress tensor for boson systems 
is symmetric. It in first order in interaction radius consists 
of two parts, namely, the term coincides with interaction 
via binary interaction potential and the term coincides with 
TPI which does not contain interaction by means binary 
potential. The first term leads to equation coinciding with 
the analogous GP equation. The second term includes ex- 
plicit expression for three particle interaction constant and 
coincides with the term analogous to the nonlinearity of 
fifth degrees in NLSE. 

In dilute alkali gases, when interaction may be consid- 
ered like scattering process, the constant of interaction 
in first order on interaction radius (FOIR) has the form 
T = —4:TTh^a/m; a- is scattering length. In TOIR ap- 
proximation there arise the second constant of interaction 
T2. In general case, parameter T2 is independent from T 
and has to considered like supplementary. In (t^ an ap- 
proximate estimate of T2 via T is considered. 

In the paper 13^ were shown the account of interac- 
tion up to TOIR leads to finding of new physical effects in 
BEC. In | 36] were found new type of solitons in BEC. Fre- 
quency dependence of eigenwave in boson-fermion mix- 
ture in TOIR is derived i n Ivl] . Therefore, generalization 
of Bogoliubov spectrum ll34ll . Issll was obtained in TOIR. 
Also, analytical dependence U!(k) for degenerate fermions 
is obtained. Usin g sy stem of QHD equations (derived in 
|7]) in works 11281 13311 was investigated dynamics of non- 
linear wave in TOIR. In ll28ll analytical solution for bright 
soliton is obtained in uniform BEC within the TOIR ap- 
proximation. In ll33ll nonlinear eigenfrequency shift in 
uniform BEC is analytically investigated within the TOIR 
approximation. 

At experimental investigation of BEC in magnetic traps 
work with number of particles N is order 10''- 10^. At 
quantum-mechanical description of such system it is neces- 
sary to solve Schrodinger equation determining wave func- 
tion that depends from 3A^ coordinate and time. Wave pro- 
cess, process of transfer, exchange by energy and momen- 
tum at interaction take place in three-dimensional physi- 
cal space. In this connection it is necessary to convert 
Schrodinger equation to equation is determined dynamic 
of functions in three-dimensional physical space. This task 
is solved with method of quantum hydrodynamics, further 
development of this method is made in this work. 

Moreover, in this work problem of finding of method al- 
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lowed to build NLSE on the base system of QHD equations 
is solved. 

Contribution of noncondensate particles was shown to be 
arise in QHD equations for BEC and degenerate fermions 

a. 

Our paper is organized as follows. In Sect.2 we derive 
QHE's for BEC from many-particle Schrodinger equation 
with two- and three-particle interaction. In Sect.3 we cal- 
culate quantum stress tensor due to three-particle interac- 
tion. In Sect.4 we derive system of QHD equations for sys- 
tem of bosons at non-zero temperature. We consider sepa- 
rate dynamic of two type of bosons, it is particles in BEC 
state and noncondensate bosons. We obtain the continu- 
ity equations and Eiler equations for each type of bosons. 
This equations include interaction between different type 
of bosons. The two- and three-particle interaction are in- 
cluded too. In Sect.5 from QHD equations finding in sect.3 
we obtain NLSE with nonlinearity of third and fifth de- 
grees. We give special attention for interaction BEC with 
temperature excited particles, but part of this presented in 
the Appendix 3. In Sect.6 we construct equations for BEC 
including two-particle interaction in FOIR, TOIR and con- 
tribution of first term of three-particle interaction. On this 
background we obtain frequency dependence of elemen- 
tary excitation on wave vector for BEC with two- and three- 
particle interaction at zero temperature limit. In Sect.7 we 
study the dispersion of elementary excitations at non-zero 
temperature. We obtain the velocity of the first and second 
sounds due to TPI. In Sect.8 brief summary of obtained 
results are presented. 



II. DERIVATION OF QUANTUM HYDRODYNAMIC 
EQUATIONS 

In the case, when we can present interaction between 
three particles Vsir.i, Fj, r^,) in the form 



(1) 



where = — r^, Uij 
teraction potential, Uijk = U{\ 



U{\ I) -is the binary in- 
TiA; I) I Tfcj |)-is 



the three particle interaction potential that does not contain 
combination of binary potentials Uij . Collisions of three 
particles in dilute gases is rare. But in the case of liquid 
helium there is interaction of more than two particles at the 
same moment. Moreover, this interaction has not form of 
scattering. 

Let us consider the system of N Bose-particles with a 
short-range potential. The Hamiltonian of the system under 
consideration has the form: 



Uijki (2) 

where = — z/iV,-is the momentum operator of the i-th 
particle, m^-is the mass of the i-th particle. Let us consider 
the three-particle interaction nonequivalent to combination 
of binary interaction consider by Uij . Interaction of three 
or more particles at the same time by means of binary po- 
tential is described by Uij and have no connection with 
TPI. 

Atoms in the traps are kept by means interaction of 
their magnetic moments with trapping magnetic field. In 
that theory this interaction implicitly take into account by 
means Vi^ext- In the magnetic field of third particle may ex- 
ist effect like Feshbach resonance |26]. Due to this effect 
scattering length of two-particles interaction is changed. It 
may be considered like mechanism of three particle inter- 
action, which is not provided by binary potentials, only ([T]l- 

The concentration of particles in the vicinity of the 
point r of the physical space is determined as the opera- 
tor (5(r — r^) averaged over the quantum-mechanical 
states: 

n(r,t)= I dRY,5{r-r,)i^+{R,t)^{R,t), (3) 

where dR = Hili dVi. 

Differentiating this function over time and using the 
Schrodinger equation with the Hamiltonian dJll, we derive 
the continuity equation, in which the current density vector 
appears in the form: 

nr,t) = JdRY.S{r-r.)^^ 



X {pf^P)+iR,t)^^JiR,t)+^P+iR,t)ip:i^)iR,t)]. (4) 



H 



1 

2m 



The momentum balance equation for the system of par- 
ticles under consideration is obtained similarly, i.e., via 
differentiation of current density (ffl and applying the 
Schrodinger equation ifsHl . i32ll . Q]. As a result, we 
obtain: 

mdtr{r,t) + dpU"^{r,t) 
= -J dr'(V"?7(r,r'))n2(r,r',t) 

- Jdr' J (ir"(9"[/(r,r',r"))n3(r,r',r",t) 

-n(r,t)V"K.t(r,t). (5) 

In the momentum balance equation, 11"^ (r, t) is the quan- 
tum tensor of the density of the momentum flux. This ten- 
sor has the form 

n"^(r,t) = J di?5]<5(r-r,) 



4 



4m, 



+ {p':i^)+iR,t){p^i;){R,t)+c.c.]. (6) 



The interaction between the particles in ^ is expressed 
through the two-particle probability density n2(r, r',t) 
normalized over N{N — 1) and having the form 

n2(r,r',t) = j dR^, <5(r - r,) 

x5{r' -r,)i;+{R,t)i){R,t), (7) 

and three-particle probability density n3(r, r', r", t), nor- 
malized over N{N — 1)(A^ — 2) and having the form: 

n3(r,r',r",t) = j dR J2 " r^) 

X 6{r' - r,)6{r" - r,)7^+(ii, t)^iR, t). (8) 

Let us represent the first term in the second member of 
equation ((Sj, i.e., the density of the interaction force of the 
particles, in the form 

-IjdRT. (5(r-r.)-^(r-r,)) 



X (Vf[/(r,,))V^+(/?,t)V'(i?,t), 



(9) 



which is possible by virtue of symmetry (antisymmetry) of 
the wave function, and let us proceed in ^ to variables of 
the center of gravity and variables of the relative distance 
of the particles: 

1 / 



(10) 



Rewrite the three-particle interaction potential in the form: 

= -d'^U{rij,rjk,rui) - d^U{r,j,rjk,rk,). (11) 

Second term on the right-hand side eq. ^ is represented 
in the form 

+5{r -Yj)d"U {rij , r^-fe , r^, ) 

+5{r - rfc)a,"C/(r„ , r,fc, r,,) j i^+{R, t)i^{R, t). 
Using (fTTl) we obtain 



+ [5{v-v,)-5{v-Tu)]d^U,,k)i^+{R,t)i;{R,t). (12) 

For the case of three particles we can use variables of center 
of mass and relative motion too. 



^ijk — + + T/c), Tij — Ti — r j , r^fc 



life 



(13) 



Since the interaction forces between the particles rapidly 
descend at distances of the order of the interaction radius, 
small \rf^\ give the main contribution to the integral (|9|. 
Therefore, in expressions ^ and (fTzl) we can replace the 
multipliers at the interaction potential by their expansion 
in series by r . We have concluded that the density of the 
interaction force for bosons with a short-range interaction 
potential can be represented in the form of divergence of 
the tensor field 5^0""^ (r,t). Here, cr"'^(r, t) is the quan- 
tum stress tensor due to inter-particle interaction. 

Therefore, the momentum balance equation will take the 
form 

dtr{r,t) + -dp{Il^^{r,t) + cT"^(r,t)) 
m 



1 



m 



n(r,t)V"K.t(r). 



(14) 



Let us now consider the tensor 11"'^ (r, t) and isolate in 
it the contributions to the momentum flow density for the 
convective and thermal motions and the purely quantum 
part. For this purpose, let us introduce velocities by formu- 
las 

y,[R,t) = —V,S{R,t), (15) 

where Vi(ii, t) is the velocity of the i-th particle, while the 
S{R, t)-phase of the wave function is 

ij{R,t) = a{R,t)expC-^^^). 

h 

Velocity field v(r, i) is determined by the formula 

j(r,t) = n(r,t)v(r,t). (16) 

Then Ui(r, R, t) = v,;(i?, t) — v(r, t) is the quantum ana- 
log of the velocity of thermal motion. Isolating the explic- 
itly thermal motion of the particles with velocities and 
the motion with the velocity v(r, t) in continuity and mo- 
mentum balance equations ( fT4l) . we come to the following 
equations: 

dtn{r,t)+V{n{r,t)y{r,t)) = 0, (17) 

mn{r, t){dt + yV)v"{r, t) + n(r, t)V''V,,tir, t) 

+ 9^(K'^(r,t) + cj"'^(r,t)+r"^(r,t)) = 0. (18) 
In equation ( fTSl) 

N 

p"/3(r,t)= dRY,Hr-r,)a\R,t)m,ufu^. (19) 
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This tensor tends to zero along with equality to zero of ve- 
locities of thermal motion Uj of the particles. Therefore, it 
has the meaning of the kinetic pressure. 

The tensor T"'^(r, t) is proportional to H"^ and has a 
purely quantum origin. For the system of numerous nonin- 
teracting particles, this tensor is 



4m 



d"d^n{r,t) 



n{r, t) 



{dMr,t)){d^n{r,t))]. 



(20) 



Therefore, in this section we have obtained general form 
of QHD equations for system of particles with short-range 
interaction. 



in. CALCULATION OF QUANTUM STRESS TENSOR 

One have the aim to calculate quantum stress tensor 
cr"'^(r, t) it is necessary to write explicit form of a°'^{r, t) 
through wave functions. 

The first terms of series in quantum stress tensor 
a°''^{r,t) gives the main contribution. Writing this terms 
alone, we have 

a"P{r,t) = -\j dR 5(r-R,,) 



I Ty- I d I I 



Where 



X J(r-R„fe)V'+(i?,t)V'(i?,t). 



^{R,t) = 'lp{...,Rijk,...,R^jk,■■■,t) 



(21) 



in the term describing binary interaction, similarly in the 
term describing three particle interaction: 

il^iR, t) = Rijk, Rijk, Rijk, t). 

Using symmetry properties both wave function of bosons 
and potential of interaction between particles we can trans- 
form formula (1211 to the view: 



C7"^(r,f) = -irr(n2(r,r',t)) j dr 



r^'r^ dU{r) 
r dr 



Tr{m{r,r',r",t))x 



ri3 I) j +r^i-i^{U{ri2,r2^, \ ri2 - rig 



(22) 



where 

rr/(r,r') = /(r,r), 

Tr/(r,r',r") = /(r,r,r). 

It is evident that the first term in the second series of 
o""'^ (r, t) (related with TPI) for the case of fermions equals 
to zero. For investigation the influence of TPI on dynamics 
of fermions one needs to use the next term of this series. 

I 



I 

Using knowledge of decomposition formulas of wave 
function and relations of orthogonality, which are pre- 
sented in appendix I, we obtain following expression for 
two-particle concentration 



n2(r, r', t) = n(r, t)n{r\ t) + |/j(r, r', t)p + p(r, r', t), 

(23) 

and for three-particle concentration 



n3(r, r', r", t) = n(r, t)n{r' ,t)n{r" , t) + n(r", t) \ p{r, r', t) \^ +n(r, t) \ p{r\ r", t) 1^ 
+n(r',t) I p{ry,t) p +[p(r, r', t)p(r', r", t)p(r", r, t) + c.c] 
+n(r', t)p{r, r", t) + n(r", t)p(r, r', t) + n(r, t)p{r' , r", t) 



+ W\ r', t) V n,{n, - 1) (^,(r, t)^,{r" , t)^*(r, t)M , t) + c.c 
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(24) 



In (|23] |. (|24] | there is not distinction functions describing 
behavior of BEC and of particles in exited states. 
Here 

^(j": *) = XI ^9<^s(''' O'/'slr, (25) 



P(r,r',t) = Xl"ff'^s(''>*)¥'9(r',i), (26) where 

9 

P(r,r',t) = - l)|^,(r, i)n(^g(r', t)!^, (27) 



where v'g(r) ^) ^re the arbitrary single-particle wave func- 
tions. 

The last two terms in formula ( l23l) represent part of stress 



6n^(r,t) + 5n(r,t)p(r,t) 



9 

p(r,t) = pB{r,t) + p„(r,t) 



and 



Trp{v,r\t) 



r,t). 



(28) 



(29) 



tensor arose via exchange interaction. Substituting expres- ^ , r- n • • i • i ■ i i i ■v- , 

ioTi • jooi r. ° ^ ^ r- 1 In the ( I28D rollowing integrals is designated through I and 

sion i23i in (1221) tor quantum stress tensor or bosons sys- — o o o o 



tem, taking into account Trp{r, r', t) = n(r, t), one ob- 
tains formula 

a"^(r,t) = -iTr^(2n2(r,t) + p(r,t)) 



Q/3 



,a[/(r) 
9r ' 



(30) 



X 



a/3 



6 J V V '^l ^2 Ti - r2 



X [/ (ri, r2, W rl + rl + 2rir2 cos il) , 



(31) 



where is angle between ri and r2 and di, 82, 83 are 
derivatives of function U on its arguments. We can see that 



Assuming that the potential satisfies the condition that 
the quantity r^U{r) tends to zero at r tending to zero and 
infinity, for T from (l30l l. by integration by parts, we obtain 



T 



drU{r), 



which coincides with the result for the interaction constant 
g found by Gross and Pitaevskii allowing for the sign in 

Taking into account the notions given after formula (l62l l. 
we obtain the following relation for the quantum tensor of 
stress of the system of bosons close to the BEC state: 



X 2ns (r, t)n^{r, t) + 2nl{r, t) + p(r, t) 



'X^ 18nB(r, tX(r, t) + 5n„(r, t)p(r, t) 



+ 6n^(r, t) + 5nB(r, t)p„(r, t) + m(r, t)j , (32) 
where 

9 

Here, we used the notations nsi^, t) for the concentration 
of particles situating in the BEC state and n„(r, t) for the 
concentration of excited particles. Notation go designates 
ground state of system of particles corresponding to BEC. 

The stress tensor that depends on inter-particle interac- 
tions contains the single-particle functions ipg{r,t). By 
these functions, expansion of the unknown N-particles 
wave function is actually performed. If we neglect the 
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inter-particle interaction, the N-particle problem will be 
reduced to the single -particle one, and the set of func- 
tions (r, t) will be determined by the single-particle 
Schrodinger equation. In this case, the momentum balance 
equation will also not contain interactions and, along with 
the continuity equation, will determine the single-particle 
wave function in the polar form. Such single -particle func- 
tions, which are simultaneously solutions of equations of 
quantum hydrodynamics of the system of noninteracting 
particles and the single-particle Schrodinger equation, can 
be used as the first approximation at calculating of the 
stress tensor. Thus obtained quantum balance equations 
will determine the set of field functions that determine the 
state of the system in hydrodynamics and, consequently, 
the single -particle wave function of the system of interact- 
ing particles. Such a wave function of spatial coordinates 
and time will determine the effective Schrodinger equation, 
which is nonlinear and integro-differential in the general 
case. The solution of such equation can be used as the sec- 
ond iteration at calculating the stress tensor, etc. 

For particles to be found in state of Bose condensation 
we can get relation for ps and ffiB 

pB{r,t) = ^ng{ng - l)|(/5<,(r, 

g 

g 

= {N\^g,{r,t)\''f = nl{r,t), (33) 

similarly: 

fnB{r,t) = ^ng{ng - l){ng - 2)x 

9 

X |v?,(r,f)|« =ni(r,t). (34) 

Using results ( l33l) , (l34l l. we obtain expression for 
(T2^„(r, t) in the form: 

<„(r,t) = -iT5"^(^2ns(r,t)n„(r,t) 

+n|(r, t) + (r, t) + p„(r, t)^ - ^x"^ (^6nl{r, t) 
+18nBir,t)nl{r,t) + 5n„(r, t)n^(r, t) 

+ buB (r, t) p„ (r, t) + n|(r , t) + m„^ . (35) 

Without exited particles quantum stress tensor takes the 
form: 

<(r,i) = -^T5"^n|(r,t) - ^x"'n|(r,t). (36) 



The corresponding momentum balance equation of the 
quantum hydrodynamics for BEC alone takes the form: 

mn(r, t){dtv''{r, t) + w'^(r, t)V^v''{r, t)) + d^p^^ir, t) 



-7^"- r, t)do,—^j==- - Tn r, t)9"n(r, t) 

-2x"^n\r,t)d^n{r,t) = -n(r, t)V"K.t(r, t). (37) 

Therefore in this section we have obtained the QHD equa- 
tions for BEC with two- and three -particle interaction. 

IV. EILER EQUATIONS FOR A BOSON SYSTEM WITH 
NONZERO TEMPERATURE 

If we interesting in separate dynamics of BEC and tem- 
perature excited bosons we can divide the concentration 
n(r, t) and current j(r, t) on two part 

n(r, t) = UBir, t) + n„(r, t) 

and 

j(r,t) =js(r,t)+j„(r,t). 

Here we consider a case there are no exchange of parti- 
cles between BEC and noncondensate component. In this 
situation we can write the continuity equations for each 
kinds of particles, for the BEC 

dtHBir, t) + V(nB(r, t)vs(r, t)) = (38) 

and, for the noncondensate bosons 

<9tn„(r, t) + V(n„(r, t)v„(r, t)) = 0. (39) 

In previous section we derive the momentum balance equa- 
tion for whole system of bosons, i.e. for the mixture of 
BEC and noncondensate particles. Now we need to obtain 
momentum balance equations for each kinds of bosons. 
For separation of contributions of different sorts of bosons 
we must to consider of derivation of quantum stress tensor 
( |23]) and ( [24] | or (|35] |, in more detailed way. 

For understanding the detail of evolution of (r, t) and 
j„(r, t) we need to consider formulas i62\ and ( l64l) in de- 
tail. The first multiplier in this formula, which has argu- 
ment (r, t) is related to the particle whose motion we con- 
sider. Another one particle wave functions are related to 
the particles that influence on dynamic of considered cur- 
rent. This is give us ability to obtain the separate equa- 
tion of dynamic atoms in BEC state and noncondensate 
ones. We consider each term in ( [62l ) and ( |64l l. If one- 
particle wave function with argument (r, t) describe BEC 
state (has subindex "B"), we put this term in momentum 
balance equation for BEC. In the case one-particle wave 
function with argument (r, t) describe noncondensed state 
we put this term in momentum balance equation for non- 
condesate particles. In this way we get following expres- 
sions for quantum stress tensors of atoms in both BEC and 
noncondensed states. 
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a 



Pb) 



and 



-X^^i^riBnl + AnnpB + nBPn + mB) (40) 



= --T5"^{nBn„ + 2nl + pr, 



--X"'^{^ril + 12nBnl+5nnPn+4:nBPn+nnPB+rh„). 

(41) 

In appendix 2 we calculate p„ and m„ and receive 

p„ = 0, m„ = 0. 

From formulas (l33l l and (l34l l we see (f) = (J*' ^) 
and iriB{T,t) = n'^(r, t). In this case we have ag^ and 
o""'^ has more simple form and expressed in terms of 
and n„. 

Obtained expression (|28]l for o-"'^(r,t) in the approxi- 
mation under consideration should be substituted into mo- 
mentum balance equation ([Tsl l. 

Using this results and formulas (|33] |. (|34] | we can rewrite 
equations (l40l l and (1411 in the form 



A^nB{r,t) 

-—nB{r,t)dc. 
2m VMM) 



+ Ix'^^d^ (gubTiI + 4n„n|^c + n^^^:) (42) 



and 



mn„ (r, t) (r, t) + i;^ (r, t) (r, t)) 



V. DERIVATION OF NLS EQUATION FOR ZERO 
TEMPERATURE 

There is used the NLSE along with QHD equations in the 
literature. From QHD equations one can obtain the NLSE 
Hill- 

Macroscopic wave function may be defined via hy- 
drodynamic variables: the concentration of particle number 
n(r,t) and potential of velocities field (p(r, t). For determi- 
nation of evolution of velocity field the integral of Cauchy- 
Lagrange following of the momentum balance equation is 
used. 

Therefore, for eddy-free motion v(r, t) = V0(r, t) and 
barotropicity condition in terms of tensor field 



mn(r, t) 



(44) 



where ^"''(r, t)-tensor of chemical potential, momentum 
balance equation has the first tensor integral: 



m m m 



, (constr^. (45) 

Continuity equation ([TTl) and momentum balance equa- 
tion (l37l l can be associated with the equivalent one-particle 
Schrodinger equation for some effective wave function 
<I>(r, t). Represent this function in the form: 



$(r, t) = ^n(r, t) exp (^^"^(/.(r, t) j . (46) 

Differentiating it by time and using equations ( [TTl ). ( |45l ) we 
take equivalent one-particle Schrodinger equation: 



ih6''^dMr,t) 



2m 



2m Vn,,{r,t) 



= -n„(r,t)V"K.t(r,t) + -T5"(nBn„ + 2nl) 

+ \xr^d^ + yiriBul + n„n|sc) ■ (43) 

where evident form of interaction constants T and x"'^ are 
presented by formulas (|30] | and (1311 correspondingly. 



- T(5"^ I $(r, t) Y -x"" I $(r, t) ^'(r, i), (47) 

under conditions x"^ = X*^"^- P"'^(r) = p(r, 
isotropic kinetic pressure and barotropicity condition: 



Vp(r,t) 



V/u(r,t), 



mn(r, t) 

where ^(r, i)-chemical potential, we obtain: 



(48) 



iHdt^{r,t) 



2m 
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+ K.t(r,t) - T I $(r,t) p -X I Hr,t) \'jHr,t). 

(49) 

That is well known Gross-Pitaevskii equation ll35ll38l - l40ll . 
with the nonlinearity of fifth degree lilljl, which arise due 
to TPI. 

The isotropic constant of three -particle interaction is: 



^ j dridr2(J^idr,+ | ri - ra | 



X 



X U{ri,r2, ^Jrf + rj + ^r^r^ cos 0)J . (50) 

Wave function $(r, t) is normalized with condition: 

y dr$(r,t)*$(r,t) = TV, 

where is number of particles in the system. 

Within this approximation the equation of momentum 
balance has form: 

mdtv'^ir, t) + ^md"v^{r, t) + md"fi{r, t) 



-x9"n2(r,t) = -5"K.t(r,t). 



(51) 



If system of particles is dense (e.g. Bose liquid) then two 
more particles can interaction simultaneously. Basic con- 
tribution in interaction three and more particles arise from 
binary interactions Uij. Part of interactions in macroscopic 
system of particles which is described via binary potential 
is the origin of first series in (|2T] | |7] and, consequently, 
this part is also the origin of the term which is proportional 
Tin (EH). 



VI. LINEAR WAVE DISPERSION IN BEC WITH TPI 
ZERO TEMPERATURE LIMIT 

In this section we have the aim both to consider disper- 
sion properties of linear wave and to compare two correc- 
tions to GP approximation. For this aim we write system 
of quantum hydrodynamic equations in third order to the 
interaction radius, for only binary interaction which was 
given in the work [7|]. Taking in account binary and three- 
particle interaction and including results of the previous 
section we obtain equation: 

mn(r, t)dtV°'{r, t) + ^mn(r, t)V"t;^(r, t) + V"/u(r, t) 



2m V?^(r, t) 2 



-T2d^An%r,t) = -n(r, t)V"Kxt(r, t). (52) 



where 



vr 
30 



dr 



Following to the article 14111 we redefined T2, in compare 
with the paper fv*]. 

Let us consider the eigenmodes, which can propagate in 
one dimensional geometry, based on the Eqs. of continuity 
( fTTl) and momentum balance (l52l) . 

We consider the small perturbation of equilibrium state 
like 



n 



Uq + Sn, v° 



+ 



(53) 



Substituting these relations into system of equations ([TTl l 
and ( l53l) and neglecting nonlinear terms, we obtain a sys- 
tem of linear homogeneous equations in partial derivatives 
with constant coefficients. Passing to the following repre- 
sentation for small perturbations 6f 



6f = f{uj,k)exp{—iLLi + ikr) 



(54) 



yields the homogeneous system of algebraic equations. 
The magnitude of concentration of BEC is assumed to have 
a nonzero value. Expressing all the quantities entering the 
system of equations in terms of the concentration of BEC, 
we come to the dispersion equation for elementary excita- 
tions 



OJ 



4m ^ 



m 



'^riQ 2x 2 1 ,2 

H riQ k^. (55) 



m 



m 



In the absence of T2 the result (l55l l is in accordance 
with the real part of solution uj{k) obtained in [i25i1 . If 
we compare terms which proportional to x and T2 we can 
see functional dependence cj(A;) for the cases three -particle 
interaction and binary interaction in third order on interac- 
tion radius. When we calculate binary interaction in third 
order to interaction radius there arise linear dependence on 
Uq. Term proportional to in (|55l l. las in GP approxi- 
mation, is linear on concentration Uq. There is additional 
dependence on in dSSl l at the expense of three-particle 
interaction. 



VII. LINEAR WAVE DISPERSION IN BEC WITH TPI 
NONZERO TEMPERATURE LIMIT 

In this section we consider a Bose particle system con- 
tained particles being in the BEC state and noncondensate 
atoms. We interested the dispersion of small amplitude el- 
ementary excitation in such system. For this aim we use 
system of equations ( [39l ). ( |38l l. ( |42] | and ( |43l l. Equilibrium 
state is described by constant values of concentrations of 
each sorts of bosons Uqb and no„. The flows in equilibrium 
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state are absent = v„ = 0. For small perturbations of 
equilibrium state we have 

n, = no. + (5n„ vf = + vf, (56) 

where i-mark B or n. Using formula (l54l) for small pertur- 
bations we get the homogeneous system of algebraic equa- 

I 



tions. A condition of existence on nontrivial solution lead 
to dispersion equation. A solution of dispersion equation 
has form 



1 fl'^k'^ 16 
±^ [^-^Kb + 9nL - 2nosno„) + y (/c"/c''x"'')'(noB + 36n^„ + 124n^5n^„ + 2AQn^Bnl^ - Snj^^non) 

4 \'/' 
+ -Tek'^k^x'^^nls + 18n;^„ - 5ngsno„ + 42nosngJ . (57) 



In the case x°'^ = X<^"^ we have k°'k^x°'^ = ^^X ^^d 
we can make k^ out of the square root. Then, all terms, 
apart from first term, are proportional to A;^. A coefficient 
at k^ it is a square of velocity of sound. Due to two sign 
at the square root we have two sound velocity, for the low 
temperature first and second sounds. 

The first sound it is a usual sound, one take place in clas- 
sic gases and in Bose system far of BEC condition. For the 
low temperature quantum velocity of first sound coincide 
sign "-" in the equation (ISTl l. Sign "+" in equation (l57l l co- 
incide to the second sound it is exist in superfiuid systems 
due to BEC. Moreover, solution with sign "+" in equation 
(ISTl l it is the Bogoliubov's mode in BEC. In comparison 



with Bogoliubov's, our result account three particle inter- 
action and influence of noncondensate particles. 

In the absence of TPI k^k^x"^ = 0' for repulsive two- 
particle interaction T < 0, from (l57l l we have 

2 h^k"" \T\k^ f 

to = - — - H X 3noB + 5no„ 

Am'^ 4m \ 



i V '^OB + 9"-on - 2noBnonj (58) 

where the quantity in bracket is positive. 
For the case T = and /c"/c^x"^ < we obtain 



9 1 1 r\j r\j\i\irXiX\l/') 9 \ 



± \ nlg + 36n4„ + I24.nlgnl^ + 240noBn3„ - Sn^^no,, (59) 



where the quantity in bracket is positive. 

VIII. CONCLUSION 

In this article we gave derivation of the NLSE for BEC 
with cubic and quintic nonlinearities from microscopic 
quantum theory. For the derivation we used method of 
quantum hydrodynamics. In original Schrodinger equa- 
tion we took into account two-particle and three-particle 
interactions. In the article, arising force field of whole 
accounted interaction in the form of divergence of ten- 



sor field, is demonstrated. We obtained the expression of 
quantum stress tensor by means of microscopic wave func- 
tion. Using its representation we derived equation of state 
for boson systems in BEC state including two- and three- 
particle interaction. Therefore, dependence of quantum 
stress tensor a°'^ (r, t) on concentration n(r, t) are demon- 
strated. Thus, momentum balance equation obtained coin- 
cides with the analogous NLSE containing nonlinearities 
of third and fifth degrees. Method of finding of NLSE 
for the wave function in the medium for the case three- 
particle interaction are developed. We showed that ob- 
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tained NLSE coincides with well-known NLSE with cu- 
bic and quintic nonlinearities. Explicit form for the con- 
stant of three-particle interaction was obtained. In partic- 
ular, here represent derivation of GP equation from many- 
particle Schrodinger equation. Tensor form of constant of 
three-particle interaction is shown. 



From derivation of QHD equations including three parti- 
cle interaction, obtained in this article, and from derivation 
QHD equations including TOIR, obtained in f^, we can 
see terms due to TOIR and TPI additional. 



In this article, frequency dependence of elementary exci- 
tation on wave vector was calculated for the case including 
FOIR, TOIR and three particle interaction. Comparison of 
contributions from TOIR and three particle interaction is 
discussed. 



A special attention we got for influence of the tempera- 
ture on dynamic of Bose particles. We made a microscopic 
derivation of two-fluid hydrodynamic equation. We con- 
fine oneself by continuity equations and momentum bal- 
ance equations for each type of bosons, i.e. particles in 
BEC state and noncondensate particles, due to two- and 
three-particle interaction, in first order on interaction ra- 
dius. In described approximation we studied a spectrum of 
elementary excitations. We obtain the dispersion for two 
waves. One of them it is generalization of Bogoliubov's 
mode, another it is low temperature approximation for the 
usual sound. 



Where (r, t\f) = {pf(r,t) — one -particle wave function. 

The first term in formula (l62l) represents the particles 
situating in two different quantum states, while the sec- 
ond term is referred to particles in the same quantum state. 
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APPENDIX 1 

For two-particle density of probability in correspon- 
dence with the definition (|7]l we have: 

n2(r,r',t) 

= A^(A^-1) J dRM-2{ni,n2,...\r,r',RN^2,t) 

X {r,r',RN^2,t\nun2,. . .), (60) 

N 

where dRN-2 = Y\. 

For three-particle density of probability we have at the 
same way: 

n3(r,r',r",t) = iV(iV-l)(Ar_2)x 
X J dRN-3{ni,n2, • • • \r,r' ,r" , RN_3,t)x 

X (r,r',r",i?Ar„3,t|ni,n2, . . .), (61) 

N 

where dR^-s = Y\_ d^k- 

Using knowledge of decomposition formulas of wave 
function {r,r' , R]s[_2,t\ni,n2 ■ ■ ■) in formulas (|33] | and 
( l34l) we use following decomposition formulas iIbtIi 



(62) 

I 

Therefore, for the particles in the BEC state, it is sufficient 
to take into account the second term in formula ( l62l l. In 
consideration of the system of bosons with the temperature 
differing from zero, where the certain number of the parti- 



(r, r', RN_2,t\ni,n2 • • •> = XI ^\f) ^ RN-2,t\ni,. . . (rif - 1), , 

/ /',/'#/ * ^ 
x{RN_2,t\ni, . . . {nf, - 1), . . . (n/ - 1), . . .) 



E J ^tI^T^ (r, t\f) {r',t\f) {R^.2,t\n„ . . . (n, - 2), . . .). 
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cles is out of the condensate, the first summand of formula 
(l62l l gives the contribution both in the case of interaction 
of excited particles with each other and in the case of their 
interaction with the particles appearing in the BEC state. 
In this case, the second term of formula (|62] | gives the con- 



The first term in formula ( l64l) represents the particles sit- 
uating in three different quantum states. The second and 
third terms in ( |64] | represents case when two particles situ- 
ating in one quantum state and one particle situating in an- 
other quantum state. The last term represents the three par- 
ticles situating in the same quantum state. For the particles 
in the BEC state, it is sufficient to take into account the last 
term in ( l64l) . For the consideration of the system of bosons 
with the temperature different from zero, where the signif- 
icant number of the particles is out of the condensate or in 
the case of the temperature when there are no macroscopic 



tribution in the interaction both between the particles ap- 
pearing in the BEC state and between the excited particles 
appearing in the same quantum state. 

Moreover using decomposition formulas of wave func- 
tion, analogously with previous case, for three particle we 
have: 



(63) 



(64) 



number of particles in BEC state, the terms in formula (l64l l 
has the following influence. The first summand of formula 
( l64l) gives the contribution both in the case of interaction of 
excited particles with each other and in the case of interac- 
tion of two excited particles with the particles appearing in 
the BEC state. The second and the third summands repre- 
sent case which one from two states may be BEC state, but 
another is always excited state. The last term of formula 
( l64l) gives the contribution in the interaction both between 
the particles appearing in the BEC state and between the 
excited particles appearing in the same quantum state. 
Using relation of orthogonality 



I 

(r,r',r",i?Ar_3,i|ni,n2 . . .) = 

/ f'J'^f f'^fJ' * ^ ^ 

x{RN-3,t\ni,.. . {ufn - 1), . . . {uf, - 1),. . . [uf - 1),. . .) 

x(i?Ar_3,f|ni, . . . (n// - 1), . . . (n/ - 2), . . .) 

+ E E ]l m7n"'i)(N-2) *l^'> (■■■ V) ('".'If) X 

x(i?w_3,i|ni, ...{uf -2),... (uf, - 1), . . .) 

X {RN-3,t\ni, . . . (n/ - 3), . . .). 

I 



(ni, . . . {nf» - 1), . . . {uf, - 1), . . . (uf - 1), . . . |ni, . . . {rig. - 1), . . . (n,- - 1), . . . (n, - 1), . . .) = 

= s{f - q)5{r - mr - g") + - q)5{r - q")6{r - q')+ 

+5{f - q'W - qW - q") + 6(1 - q')5{f' - q")5{f" - q) + 
+6{f - q"nf - q)5{f" - q') + 5{f - q")5{f' - q')5{f" - q) 

{ni, . . . {uf, - 1), . . . (n/ - 1), . . . |ni, . . . {uq, - 1), . . . (n, -!)>•••) = 

= Kf-qW-q') + Kf-(i'W-q). 



and we obtain the expression for two- and three-particle con- 

(ni, . . . (n/ - 2), . . . |ni, . . . (n^ - 2), . . .) = 5{f - q), 
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centration presented in Sect.3. 



APPENDIX 2 



Details of calculations of p„(r, t) and m„(r, t) are pre- 
sented here. The quantity p„(r, t) and m„(r, t) reads 



P„(r,f) = rig{ng - l)\^g{Y,t) 



and 



m„(r,t) 



5]n,(ng-l)K-2)|^,(r,t)|^ 



The quantum number g indicating quantum state is the 
energy e. Therefore, n^- is the Bose distribution function, 
it's evident form is 

= = e(^.-M)V-l ' ^^^^ 
where T-is the temperature, ^u-is the chemical potential. 
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and 



2m, 



Analogously, we write for m„(r, t) 



where Qq is a label for state with lowest energy and coincide 
to the BEC state. 

We can calculate p„(r, t) and m„(r, t) approximately. 

We use method of calculation of correlations described 
in papers |7] and [31]. At the first step we suppose the 
particles is free. In this case motion of particles described 
by plane- wave. 

For the case of plane waves 

^p{r,t) = -^ex.p(^-^{ept-pr)j, (65) 

where Sp- is the energy of wave with momentum p. In 
this case we have |93g(r, t)| = 1/VV. Consequently, for 
p„(r, t), in the many-particle system we get 



V2: 



I V£den,{n, - 1) 



where g is the Lande factor. 



: / yeaen^. 



(66) 



gm 



3/2 



J y/eden,{ne - l){n^ - 2) 



gm 



3/2 



(68) 



The quantity p„ (m„) expressed in terms of sum of 
squares (cubes) of occupation numbers of quantum states 
(l66l) . (l68l l. Since sum of squares (cubes) is much less than 
square (cube) of sum, we obtain p„ <^ n^, m„ <C n^. 
Thus, we neglect by p„ and m„ in equation of quantum 
hydrodynamics (1411 . In formal way, we can write 

p„ = 0, m„ = 0. 



APPENDIX 3 

Analogously to section (5) we can enter the macroscopic 
one-particle wave function for each subsystem of bosons, 
for BEC 



^B{r,t) = ^JnBir,t)exp\^-m(t)B{r,t) ] (69) 
and for noncondensate bosons 



^>„(r,t) = ^n„(r,t)exp(^-m(/>„(r,t)j. (70) 

Using equations ([381), ^9^, (|42l) and (|43l) we have gotten 
NLSEs for each type of bosons. For particles in BEC state 
the NLSE arise in the form 



/ ^^V^ / 1 1 [ Tin 

ihdt^B{r,t) = — h lJ-B{r,t) + Ve^t{r,t) - T + -n„ + - / —dni 

\ 2m \ 2 2 J Ub 
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x(^nl + 4nl + 4: j {^dne + ^ngdn^ + ^n^dnB)j^<^Bir,t) (71) 



and, for noncondensate component NLSE reads 



drij- 
2n„ 



ihdt^nir, t) = I — — + /i„(r, t) + Kxt(r, t) - T I -n^ + 2n„ + / ■ 

- X ( + 8nBn„ + + ^ y (^!^£±i^^ risdn,, j J $„(r, t), (72) 



where 

riB =1 $s(r,t) 1^ 

and 

n„ =1 $„(r,t) |2 . 

Wave functions <I>B(r, t) and <^„(r, t) are normalized 
with conditions: 

J dr<^>Bir,tr<^>B{r,t) = Nb, 

J dr$„(r,f)*$„(r,t) = A^„. 

where Nb and A^„ are number of particles in the BEC con- 
dition and in the noncondensate particles. 
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